Band structure asymmetry of bilayer graphene revealed by infrared spectroscopy 
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We report on infrared spectroscopy of bilayer graphene integrated in gated structures. We ob- 
served a significant asymmetry in the optical conductivity upon electrostatic doping of electrons 
and holes. We show that this finding arises from a marked asymmetry between the valence and con- 
duction bands, which is mainly due to the inequivalence of the two sublattices within the graphene 
layer. From the conductivity data, the energy difference of the two sublattices is determined. 



Recently there has been unprecedented interest in 
carbon-based materials due to the discovery of graphene 
[l|. Among all carbon systems, bilayer graphene stands 
out due to its remarkable properties such as the forma- 
tion of a tunable band gap between the valence and con- 
duction bands [3t3|[j|[5|: a property not attainable in 
common semiconductors. The vast majority of previous 
experimental and theoretical studies of bilayer graphene 
assumed a symmetric band structure that is governed 
by the interlayer coupling energy 71 . This is in contrast 
with a significant electron-hole asymmetry observed in 
cyclotron resonance ^j and cyclotron mass experiments 
[^. Several theoretical proposals have been put forward 
to explain these results 2][2l]- Yet, the microscopic origin 
of the observed effects remains unknown. 

Here we present the first investigation of the opti- 
cal conductivity of bilayer graphene via infrared spec- 
troscopy. We observed dramatic differences in the evolu- 
tion of the conductivity for electron and hole polarities 
of the gate voltage. We show that small band parame- 
ters other than 71 give rise to an asymmetry between the 
valence and conduction bands, in contrast to the com- 
monly assumed symmetric band structure. The system- 
atic character of our IR data enables us to extract an 
energy difference between the A and B sublattices within 
the same graphene layer (Fig 1(b)) of Sab «18meV. We 
analyze some of the implications of these findings for 
other properties of bilayer graphene. 

Infrared (IR) reflectance R(a;) and transmission T{uj) 
measurements were performed on bilayer graphene sam- 
ples on Si02/Si substrate [6|] as a function of gate volt- 
age Vg at 45K employing synchrotron radiation, as de- 
scribed in iS]. We find that both R(a;) [i,[io| and T(w) 
spectra of the bilayer graphene device can be strongly 
modified by a gate voltage. Figure 1 shows the trans- 
mission ratio data at several voltages normalized by data 
at the charge neutrality voltage Vcn- T(V)/T(VcAr), 
where Ycn is the voltage corresponding to the minimum 
DC conductivity, and V= Yg-YcN- The T(V)/T(Vcn) 



spectra are dominated by a dip at around 3000 cm~^, 
the magnitude of which increases systematically with 
voltage. Apart from the main dip, a peak was ob- 
served in the T(V)/T(VcAr) data below 2500 cm~i, 
which shifts systematically with voltage. This latter fea- 
ture is similar to the T(V)/T(VcAr) data for single layer 
graphene [8|. The gate-induced enhancement in trans- 
mission (T(V)/T(VcAr)>l) below 2500 cm^i and above 
3500 cm^^ implies a decrease of the absorption with volt- 
age in these frequency ranges. 

The most informative quantity for exploring the quasi- 
particle dynamics in bilayer graphene is the two dimen- 



sional (2D) optical conductivity ai (w) -I- ia2 (to) [8l|[ll|. 
First, we extracted the optical conductivity at Ycn from 
the reflectance data (not shown) employing a multilayer 
analysis of the device [8|[ll|. We find that txi (w, Vcn) 
has a value of 2 * {ne^ /2h) at high energies, with a pro- 
nounced peak at 3250 cm~^ (inset of Fig 2(b)). This 
observation is in agreement with theoretical analysis on 
undoped bilayer graphene[12|[13J|14l|. Our high energy 
data agree with recent experiments in the visible region 
[l5(. The peak around 3250 cm~^ can be assigned to the 
interband transition in undoped bilayer near the inter- 
layer coupling energy 71. 

An applied gate voltage shifts the Fermi energy E^ 
to finite values leading to significant modifications of the 
optical conductivity. The ci (w, V) and (J2 (1^, V) spectra 
extracted from voltage-dependent reflectance and trans- 
mission data [8| are shown in Fig 2. At frequencies be- 
low 2500 cm~^, we observe a suppression of tJi (w, V) be- 
low 2 * {ire^ /2h) and a well-defined threshold structure, 
the energy of which systematically increases with volt- 
age. Significant conductivity was observed at frequencies 
below the threshold feature. These observations are sim- 
ilar to the data in single layer graphene [8| . The thresh- 
old feature below 2500 cm~^ can be attributed to the 
onset of interband transitions at 2Ei^, as shown by the 
arrow labeled ei in the inset of Fig 2(a) and (b). The 
observed residual conductivity below 2Ei? is in contrast 



to the theoretical absorption for ideal bilayer graphene 
13j 14| that shows nearly zero conductivity up to 2F,p. 
Similar to single layer graphene, the residual conductiv- 
ity may originate from disorder effects [l3| or many body 
intertactions [8||. Apart from the above similarities, the 
optical conductivity of bilayer graphene is significantly 
different from the single layer conductivity. First, the en- 
ergy range where the conductivity cti (w, V) is impacted 
by the gate voltage extends well beyond the 2Ei? thresh- 
old. Furthermore, we find a pronounced peak near 3000 
cm~^, the oscillator strength of which shows a strong 
voltage dependence. This peak originates from the in- 
terband transition between the two conduction bands or 
two valence bands (inset of Fig. 2a) 
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The voltage dependence of the Fermi energy in bilayer 
graphene can be extracted from (T2 (w, V) using a similar 
procedure as in [8|]. In order to isolate the 2Ei? feature, we 
fit the main resonance near 3000 cm~^ with Lorentzian 
oscillators and then subtracted them from the experi- 
mental 0-2 (^, V) spectra to obtain Cj ('^j ^)- The lat- 
ter spectra reveal a sharp minimum at uj=2Fif (Fig 2(c)) 
in agreement with single layer graphene [8|. Figure 3a 
depicts the experimental 2Ei;' values along with the the- 
oretical result in [5|. Assuming the Fermi velocity vp in 
bilayer graphene is similar to that in single layer graphene 
(vi?=l.lxlO^ rn/s), we find that our data can be fitted 
with 7i=450±80meV. Equally successful fits can be ob- 
tained assuming the Fermi velocity and interlayer cou- 
pling in the following parameter space: Vi;'=1.0-l.lxlO^ 
m/s and 71 =360-450 meV. Previous theoretical and ex- 
perimental studies showed that an applied gate voltage 
opens a gap A between the valence and conduction bands 
[2] ;3][4][4 Because A{V) is much smaller than 2Ef(V) 
for any applied bias in bottom-gate devices [5|, it has 
negligible effects on the experimentally observed 2Fip{V) 
behavior. 

The central result of our study is an observation of a 
pronounced asymmetry in evolution of the optical con- 
ductivity upon injection of electrons or holes in bilayer 
graphene. Specifically, the frequencies of the main peak 
i^peak in CTi (i^, V) are very distinct for Ep on the elec- 
tron and hole sides, as shown in Fig 3(b). In addition, 
^peak on the electron side shows a much stronger volt- 
age dependence compared to that on the hole side. All 
these features are evident in the raw data in Fig.l, where 
the resonance leads to a dip in T(V)/T(VcAf) spectra. 
These behaviors are reproducible in multiple gated sam- 
ples. Such an electron-hole asymmetry is beyond a sim- 
ple band structure only taking 71 into account, which 
predicts symmetric properties between electron and hole 
sides. 

We propose that the electron-hole asymmetry in our 
(Ti (w, V) data reflects an asymmetry between valence 
and conduction bands. Such an asymmetric band struc- 
ture arises from finite band parameters Sab and V4, 
where Sab (denoted as A in 16| 13|) is the energy dif- 



ference between A and B sublattices within the same 
graphene layer, and V4=74/7o. 74 and 70 are defined 
as interlayer next-nearest-neighbor coupling energy and 
in-plane nearest-neighbor coupling energy, respectively 
Iq 13j . We first illustrate the effects of Sab and V4 on 
the energy bands of bilayer graphene Ei(k) (1=1,2,3,4), 
which can be obtained from solving the tight binding 
Hamiltonian Eq (6) in Ref. [l3J. We find that finite val- 
ues of Sab and V4 break the symmetry between valence 
and conduction bands, as schematically shown in the in- 
set of Fig 2(a). Specifically, Sab induces an asymmetry 
in El and E4 bands such that Ei >-E4 at k=0, whereas 
V4 induces an electron-hole asymmetry in the slope of the 
valence and conduction bands. With finite V4, the bands 
El and E2 are closer and E3 and E4 are further apart at 
high k compared to those with zero V4 value. 

Next we examine the effects of Sab and V4 on ai (uj,V). 
It was predicted theoretically IJ] that the main peak 
in (Ti (w, V) occurs in the frequency range between 
two transitions labeled e2 and Cs as shown in the in- 
set of Fig 2(a) and (b). Here e2=-E4(k=0)-A/2 
and e3=E3(k=kF)— E4(k=ki?) for the hole side, and 
e2=Ei(k=0)-A/2 and e3=Ei(k=kF)-E2(k=kF) for the 
electron side|14|. with A defined as the gap at k=0. For 
zero values of Sab and V4, e2 and e3 transitions are iden- 
tical on the electron and hole sides. The finite values 
of Sab and V4 induce a significant inequality between e2 
and e3 on the electron and hole sides. We first focus on 
the low voltage regime, where u!peak=&2=G3- Because vp 
and V4 always enter the Hamiltonian in the form of v^k 
and V4k products [l3|, these terms give vanishing contri- 
butions at low V, where k goes to zero. Consequently, 
^peak value at low bias is solely determined by 71 and 
Sab, with Wpeafe =7i + ^ab and ujpeak=li - Sab for the 
electron and hole sides, respectively. At Ycn (OV), inter- 
band transitions between the two conduction bands and 
the two valence bands are both allowed, which leads to 
a broad peak centered between 71 -I- Sab and 71 — Sab 
(Fig 3(b)). From the two distinct low voltage ujpeak val- 
ues on the electron and hole sides shown in Fig 3(b), 
the values of 71 and Sab can be determined with great 
accuracy: 7i=404±10meV and (5y!iB=18±2meV. There- 
fore, the (Ti (w, V) data at low biases clearly indicates an 
asymmetry between valence and conduction bands in bi- 
layer graphene due to finite energy difference of A and B 
sublattices. 

In order to explore the V dependence of ojpeak and 
the width of the main peak in cti (a;,^) , Tpeak, we plot 
the e2 and e3 transition energies [lj| as a function of 
V (Fig. 3b), using the gap formula A(V) in [5| IT]] and 
our calculated asymmetric dispersion Ei(k) (1=1,2,3,4) 
[l|, with VF=l.lxlOSm/s, 7i=404meV, (5AB = 18meV, 
and for both V4 = and V4 = 0.04. We find that e2 does 
not depend on V4 [17| , whereas e3 is strongly affected by 
V4. With a finite value of V4 («0.04), an assignment of 
^peak to (e2-l-e3)/2 appears to fit our data well on both 



electron and hole sides. Nevertheless, larger separation 
of e2 and ea on the hole side is inconsistent with the rel- 
atively narrow peak in cti (uj, V) for both electron and 
hole injection with nearly identical width. Yet the finite 
value of the V4 parameter is essential to qualitatively ac- 
count for the voltage dependence of tOpeak, because with 
V4«0 iOpeak foUows 62 and ea on the electron and hole 
sides (Fig 3b), respectively, eluding a consistent descrip- 
tion. A quantitative understanding of the V dependence 
of LUpeak and Tpeak is lacking at this stage. Our results 
highlight a need for further experimental and theoretical 
investigation of V4 including its possible V dependence. 

We stress that 71 and 6ab are determined from the 
low bias (low ki;-) data. Therefore the values of 71 and 
Sab reported here do not suffer from the currently in- 
complete understanding of V dependence of tOpeak and 
Fpeafc discussed above. The 71 value (404±10meV) is di- 
rectly determined from measurements of transitions be- 
tween the two conduction bands or valence bands. It 
has been predicted theoretically that the band structure 
of bilayer graphene as well as the parameters 71 and Yp 
can be strongly renormalized by electron-electron inter- 
actions [7J. The 71 value inferred from our data is close 
to theoretical estimates of the renormalized 71 [7] . 

IR measurements reported here have enabled accurate 
extraction of 5ab in bilayer graphene free from ambigu- 
ities of alternative experimental methods. Interestingly, 
the energy difference between A and B sublattices Sab 
in bilayer graphene (ISmeV) is much greater than that 
in graphite {5ab ~8meV) [16[ . Such a large value of 5ab 
in bilayer is very surprising. We propose that this obser- 
vation stems from different interlayer coupling between 
the B sublattices in bilayer graphene and graphite. In bi- 
layer, the direct interlayer coupling between Ai and A2 
(Fig 1(b)) considerably enhances the energy of A sub- 
lattices due to Coulomb repulsion between the tt orbits. 
However, the sublattices Bi and B2 are not on top of each 
other as shown in Fig 1(b) and thus are more weakly 
coupled. Therefore, the energy of the B sites is lower 
than that of the A sites within the same layer, leading 
to a large Sab in bilayer graphene. On the other hand, 
in graphite the B sublattice in the third layer B3 is on 
top of Bi, and that in the fourth layer B4 is right above 
B2. The coupling between the B sites in the next nearest 
neighbor layers (Bi and B3, B2 and B4, etc) increases the 
energy of the B sites compared to that in bilayer, giving 
rise to a smaller Sab value. 

The asymmetry between valence and conduction bands 
uncovered by our study has broad implications on the 
fundamental understanding of bilayer graphene. An 
electron-hole asymmetry was observed in the cyclotron 
resonance [6| and cyclotron mass experiments [2| in bi- 
layer, both of which have eluded a complete understand- 
ing so far. Our accurate determination of finite values 
of Sab and V4 calls for explicit account of the asym- 
metric band structure in the interpretation of the cy- 



clotron data. Moreover, the different Sab values in bi- 
layer graphene and graphite reveal the importance of in- 
terlayer coupling in defining the electronic properties and 
band structure of graphitic systems. 

During the preparation of this paper, we became aware 
of another infrared study of bilayer graphene by A.B. 
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The gap formula A(V) in [5| did not take into account 
5ab and V4. However, we find that finite 5ab and V4 
values have no efi'ect on the gap A. Specifically, Sab pri- 
marily modifies the Ei and E4 bands, while leaving un- 
changed the gap A between the E2 and E3 bands at k=0, 
as shown in Fig 2(a). In addition, V4 always appears in 
a term V4k in the Hamiltonian;L|], therefore it has zero 
efi'ect on A. 
In the calculation of the energy bands Ei(k) (i=l,2,3,4), 
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we used the approximation A=0, which can be justified 
for the purpose of estimating e2 and es. The gap A is very 
small (<80 meV) in the voltage range studied in our work 
[^, and does not affect the higher energy bands Ei(k=0) 
or E4(k=0) and therefore the value of e2. Note that A=0 
is only assumed when calculating Ei(k), but not in the 
A term in the expression of e2. Moreover, E2 and E3 
bands are modified by the gap only at energies below 
A/2. Because Ef is much larger than A/2 under applied 
voltage [3], E2(k=ki7)and E3(k=ki7) are not affected by 
A. Therefore, a finite gap does not modify the value of 
ea compared to that with A=0. 



FIG. 1: (color online) T(V)/T(Vciv) spectra of bilayer 
graphene. (a) and (b): data for Ef on the hole side and elec- 
tron side. Inset of (a) : a schematic of the device and infrared 
measurements. Inset of (b): a schematic of bilayer graphene. 
The solid lines indicate bonds in the top layer AiBi, whereas 
the dashed lines indicate bonds in the bottom layer A2B2. 
The sublattice Ai is right on top of the sublattice A2. 
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FIG. 2: (color online) The optical conductivity of bilayer 
graphene. (a) and (b): a\ {oJ-,V^ data for Ef on the hole 
side and electron side, (c): ctj (a;, V) spectra in the low 
frequency range, after subtracting the Lorentzian oscillators 
describing the main resonacne around 3000 cm~^ from the 
whole (J2 (tJ, V) spectra. Inset of (a): Schematics of the band 
structure of bilayer with zero values of &ab and V4 (red) and 
finite values of &ab and V4 (black) , together with allowed in- 
terband transitions. Insets of (b): ai (oj, V) at OV (Vciv) and 
40V on the hole side with assignments of the features. 
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FIG. 3: (color online) (a) Symbols: the 2Ef values extracted 
from the optical conductivity detailed in the text. The error 
bars are estimates of the uncertainties of ctj {lo, V) spec- 
tra in Fig 2(c). Solid lines: the theoretical 2Ef values us- 
ing vf = 1.1x10^ m/s and 7i=450meV. (b) Solid symbols, 
the energy of the main peak ujp^ak in the cri (uj, V) spec- 
trum. Open symbols: the energy of the dip feature lomp 
in the T(V)/T(Vc]v) spectra. Note that ujpeak in (Ji {ijJ,V) 
is shifted from ujdip in the raw T(V)/T(Vciv) data with an 
almost constant offset, which is due to the presence of the 
substrate. Solid lines: theoretical values of the transitions at 
e2, es and (e2-l-e3)/2 with VF=l.lxlO®m/s, 7i=404meV and 
5As = 18meV and V4=0.04. Red dashed lines: es with similar 
parameters except V4=0. 



